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Abstract 

Information theory is used to perform a thermodynamic study of non 
equiUbrium anisotropic radiation. We limit our analysis to a second-order 
truncation of the moments, obtaining a distribution function which leads 
to a natural closure of the hierarchy of radiative transfer equations in 
the so-called variable Eddington factor scheme. Some Eddington factors 
appearing in the literature can be recovered as particular cases of our two- 
parameter Eddington factor. We focus our attention in the study of the 
thermodynamic properties of such systems and relate it to recent nonequi- 
librium thermodynamic theories. Finally we comment the possibility of 
introducing a nonequilibrium chemical potential for photons. 
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1 Introduction 

The study of radiation hydrodynamics has proven to be of great interest in 
astrophysics, cosmology and plasma physics. The radiative transfer equation for 
the specific radiation intensity /(r, t, u, Q) = hue n{f, t, u, Q), where n{f, t, u, Q) is 
the occupation number of photons with frequency u moving in direction fl, is in 
many practical situations too involved to be solved analytically. What is usually 
done is to consider the equations for the moments of /(r, t, z/, f2) up to a given 
order m [|, 0]. However, due to the dependence of the equation for the moment 
m on the moment m + 1, one needs to introduce a closure relation. If only the 
energy density e (m = 0) and the energy flux Je {m = 1) are considered, one 
must introduce a closure relation for the pressure tensor P^; (m = 2). Hence, in 
this approximation the relevant physical quantities are the angular moments of 
the intensity (note that in the following we consider as variables the moment of 
the photons p = Qhu/c = pc/c instead of the frequency u and the solid angle Q) 

e(f, t) 

namely the energy density, the energy flux and the pressure tensor. It is also 
convenient to define the following normalized quantities 



^ J pcn{f,p,t)d^p, 



1 



pccn{r,p, t)d^p 



^3 , pcn{f,p,t)d^p, 



(1) 

(2) 
(3) 



The closure relation is performed by the introduction of the so-called Eddington 
factor xifs), defined as the eigenvalue of the pressure tensor corresponding to 
the eigenvector n (unitary vector in the direction of the energy flux), i.e. 

TEU^Xn- (5) 
This definition leads to the relation 

= ^^I+^-^nn, (6) 

where I is the identity matrix. In the limit of isotropic radiation (Eddington 

limit) x(0) = 1/3, while in the free streaming case x(l) = 1- 

In addition, one can define in an analogous way the angular moments of the 

occupation number n{f,p, t), namely the total photon number density, the photon 

flow and the flux of the particle flow 

1 f°° 

P(^^) = ^y^ n{f,p,t)d^p, (7) 

Mr, i) = Y^l Mr, P, t)d'p, (8) 

1 r°° 

PAr(r, t) = j^l ^n{r, p, t)d'p, (9) 
and the corresponding normalized quantities 

fN-k T. = ^. (10) 
pc pc^ 

Whenever the angular and the frequency dependence of the radiation intensity 
(or the occupation number) factorize, the two sets of moments are not indepen- 
dent, but verify the relations 

Ie — Jn, T^e — Tjv. (11) 



This fact is implicitly assumed in most papers on the subject 0]. In many in- 
stances, however, the frequency is not integrated 0, |] and, thus, in this situation 
the two sets of moments also coincide. However, all the previous quantities would 
be frequency-dependent, except, again, in the case that radiation intensity fac- 
torizes into a frequency-dependent part and an angular-dependent part. 

There is a great amount of different Eddington factors in the literature in- 
troduced following physically different approaches (see for a review). Among 
the different approaches to obtain variable Eddington factors, some authors have 
used a maximum entropy principle, both from a macroscopic p, ^ and from a 
microscopic point of view 0, ^, ^. From a macroscopic viewpoint, balance equa- 
tions for the energy density and the energy flux have been considered, and an 
entropy principle has been introduced to exploit these constraints, leading to the 
so-called Lorentz's Eddington factor that Levermore had previously obtained. 
In PI the same result was obtained from the point of view of information theory 
with fixed energy flux. However, as pointed out in 0, this anisotropic Eddington 
factor corresponds to a Lorentz transformation of the equilibrium one so that 
the anisotropy is only due to the fact that the system is being observed from a 
moving reference frame. Thus, this Eddington factor does not describe a real out 
of equilibrium situation. 

On the other hand, both Minerbo and Fu Q have obtained different Ed- 
dington factors, considering as a constraint the photon flux instead of the energy 
flux. In this reference frame where the anisotropy disappears can not 
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exist. 

In this paper, we apply information theory to generahze the two situations 
mentioned above, by considering simultaneously the energy flux and the photon 
flow as independent variables. Hence, the previous results can be recovered as lim- 
iting cases, while some new situations can be analyzed. In addition, information 
theory allows a complete thermodynamic study of radiation out of equilibrium. 
The closure scheme previously described, in which the pressure tensor is written 
as a function of the energy flux, departs from the hypothesis of local equilibrium, 
which implies that the distribution in momentum space is locally (i.e., for each 
position) the same as an equilibrium distribution. However, radiation usually 
has a distribution markedly different from a black-body distribution, and this 
hypothesis must, therefore, be abandoned in radiative transfer problems [|T^. 



The structure of classical nonequilibrium thermodynamics (for example, bilinear 
forms for the entropy production rate) also presents a lack of consistency for 



radiation When the radiation field is strongly anisotropic, the mean free 

path of the photons is large and the set of macroscopic quantities describing the 
local state of the photon gas arises from interactions occurring over large regions, 
whereas when the photon gas is in equilibrium the interactions which thermalize 
photons and matter take place in a specific volume. Thus, to describe the radi- 
ation gas, the appropriate set of quantities must contain information about its 
angular distribution [Q. 

This can be done in the closure scheme described above by including the 



fluxes among the set of tliermodynamic variables witli tlie lielp of information 
theory. This formahsm can be used as an heuristic method to find a distribution 
function consistent with the information available about the system. We will be 
able to analyze also the influence of the dissipative fluxes in the nonequilibrium 
equations of state. This procedure to study nonequilibrium equations of state 



has already been used in P, [TT|, |1^ in the case of an ideal gas and in to study 
heat conduction in a boson gas. Indeed, from the point of view of nonequilibrium 
thermodynamics, the meaning of the fundamental thermodynamic quantities in 
nonequilibrium states is a basic challenge, so that it deserves attention from all 
possible points of view. 

The plan of the paper is as follows. In Section 2, we apply the information 
theoretical formalism to the study of radiation (within the low-occupation number 
approximation) under an energy and a photon flux, considered as independent 
variables. Hence, we can obtain a two-parameter Eddington factor depending on 
both fluxes and generalized equations of state. We can also recover Lorentz's and 
Minerbo's Eddington factors in the proper limits. In Section 3 we analyze the 
possibility of introducing a nonequilibrium chemical potential for photons and its 
physical consequences. The last section is devoted to review the main conclusions 
of the paper. 
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2 Anisotropic radiation under energy and par- 
ticle fluxes 

Information theory was introduced in 1957 by Jaynes in statistical mechanics 
0, ^ in order to provide a probabihstic basis to equihbrium thermodynam- 
ics. However, although the foundations of the use of the informational entropy 
functional for nonequilibrium situations are far from being trivial, it has also 



been applied to nonequilibrium situations (see [T^ and references therein). The 
method asserts that the steady state of a system, defined by the values of a set 
of macroscopic constraints, is the most (microscopically) disordered state com- 
patible with these constraints, while disorder is measured by means of Shannon's 
entropy, which is defined as follows: let M be the number of microstates compat- 
ible with the macroscopic constraints acting on the system and pi the probability 
of a given microstate i. The informational entropy is given by: 

AT 

S = -kBY.P^^^P^^ (12) 
i=l 

where /c^ is the Boltzmann's constant and providing that the normalization con- 
dition Yl,iPi = 1 is fulfilled. By maximizing ([T^) subject to the constraints (mean 
values of extensive quantities controlled in a given experiment) the probability 
of each microstate is obtained. This method realizes a probability assignment 
which is as unbiased as possible and avoids any unwarranted assumption beyond 
the information contained in the constraints. The probabilities Pi corresponding 
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to equilibrium ensembles can be easily derived from this postulate of maximum 
entropy as shown in [ 1^ . 

In nonequilibrium situations, one considers that the generalized entropy func- 
tional ( |12|) depends both on the equilibrium and nonequilibrium constraints act- 
ing on the system and the corresponding probabilities pi are obtained by the 
maximization of this entropy. These nonequilibrium constraints may be, for in- 
stance, the heat flux or the viscous traceless pressure tensor. 

If quantum systems are considered, the statistical entropy can be calculated 



in terms of the occupation number nj according to [IS 



i=l 



a) ln(l — a — 

Hi a Qi 



(13) 



where gi is the degeneracy, and a = 0, +1, —1 for classical particles, fermions and 
bosons, respectively. 

Our purpose in this paper is to generalize some previous works [H 11 H 01 by 
means of information theory and derive a more general form for the Eddington 
factor. We will study the case of a radiation gas in the low occupation number 
approximation (in order to obtain analytical expressions) submitted both to an 
energy flux an a particle flux. The reason is that, in order to obtain a non purely 
advective energy flux, as in the Lorentz case and to generalize the study performed 
by Minerbo, we must consider together the constraints of fixed energy density 
e, energy flux Je and particle flow J^. If Je and Jn are taken as independent 
variables, it is possible to demand that the particle flow be null in order to 



eliminate from the energy flux any advective contribution so that it reduces to a 
pure heat flux. Notice that the distribution function that maximizes the entropy 
can no longer be an equilibrium one, as there is no reference frame in which one 
could find an equilibrium system simultaneously at rest (i.e., with no photon 
flow) and with an energy flux. In addition, from the general expressions obtained 
in this section, one can recover in the appropriate limit, expressions previously 
obtained in the literature for radiation submitted to an energy flux and to a 
particle flux. 

Let us note that several reasons can induce a non Planckian distribution for 
a photon gas coupled to matter in a nonequilibrium state submitted to high gra- 
dients or rapidly varying fluxes. On the one hand, photons can interact weakly 
with matter in the time scale over which the flow variables change. For this 
reason, the temperature of the photon gas can be different from the local equi- 
librium temperature of the matter. In these situations, a Stefan-Boltzmann-like 
law, namely aT^, may be used with a temperature Tr different from the local 
equilibrium temperature of the matter. This situation appears for example 



in the so-called diluted radiation [0, of interest in photo voltaic devices. On 
the other hand, if we consider matter submitted to large temperature gradients 
and the length scale of its interaction with photons is large in comparison to the 
scale of variation of temperature in matter, the distribution of photons can be 
anisotropic. The previous examples are attempts to include these effects in the 
statistical distribution of photons and in the Eddington factor. 
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The distribution function to consider will then be 



2 exp —Ppc — I-pcc — K-c (14) 



exp (/3pc + / ■ pec + K ■ cj — 1 
where /, K are the Lagrange multipliers related to energy flux and particle flux 
respectively and the factor 2 is related to the two possible polarizations of pho- 
tons. If / = the constrained flux is the particle flux and one recovers the results 
obtained by Minerbo while when one takes K = 0, the constrained flux is the 
energy flux and the Lorentz limit is recovered. Any integral quantity defined in 
terms of this latter approximated distribution converges provided that (3 > \I\c. 
In addition, in order to simplify the calculations, we will assume that / and K 
are parallel. With these requirements we can find any thermodynamic quantity 
in terms of the special functions defined by: 

which verify the useful properties: 

^^i|*.-*„_J, (16) 

In particular, the Lagrange multipliers K, I can be related with the dissipative 
fluxes by using (|14D in Eqs.(|D,(|D, leading to new nonequilibrium equations of 
state: 

247rr 1 

^''=MW'b^^'^'''^^~^'^'''^^^' ^^^^ 
11 



J, 



N 



(19) 



{(3hc)^ b 

where we have defined a := \K\c and b := \I\c/P, while the caloric equation of 
state can be computed from (P and is given by: 



247r 



;^'4(a,&). 



(20) 



Note that the presence of dissipative fiuxes is seen to modify the Stefan-Bohzmann 
law. The density of photons can be also computed from (|^) and is given by: 



Stt 



{hcpy 



r^3(a,6). 



(21) 



The reduced fiuxes /e and /a? are given by: 



fE 



N 



,^4(a,&) 
^2ia,b) 



- 1 



(22) 



(23) 



_^3{a,b) 

so they do not coincide in general. The Eddington factor calculated from (|^) and 
(ID is given by: 



X 



1 

62 



1 -2 



^4(0,6) ^4(0,6) 



[/iV -fE + bfNfE] 



(24) 



and, thus, it is given in parametric form as a function of both fiuxes, /e and /at 
by Eqs. (^),(^) and (p^. The behavior of this Eddington factor for low fiux 
values is analyzed below, but let us note that in equilibrium (/at = /b = 0), the 
isotropic Eddington factor x = 1/3 is recovered. Now, we discuss Eq. 
some important particular cases. 



m 
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We start by considering the situation of a pure heat flow, which corresponds 
to energy transport without net mass flow, i.e. we impose that Jn = 0. From 



Eq. (|1^), we can easily observe that this condition is simply given by: 

^2{a,b) = ^,{a,b), (25) 



thus a and b are no longer independent variables. Using equations (|22|) and (P5|) 
one can express a and 6 as a function of Je- In this case of pure heat flux, the 
Eddington factor adopts a simpler form, 

1 



^ -^2{a,b) 



~. (26) 



^4(a,&). 

Here we recall that first the photon flux /at is constrained and then is settled to 
a posteriori. 

If the photon flux fj^ is unconstrained, which corresponds to the limit a = 0, 
/at is related to fE by: 

1 



= -(2-^/4-3/1), (27) 

lE 

and one recovers the so-called Lorentz Eddington factor, namely 

5 2 



X=3-3V4-3/|. (28) 

This Eddington factor was proven to be corresponding to an equilibrium mov- 
ing system in 0: the constraint of given energy flux Je (or equivalently, given 
momentum P, as Je = (?P) can be realized with a (local-)equilibrium moving 
system if one does not constrain the particle flow of the system. 
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In the limit 6 — > (fixed particle flux but unconstrained Je) we obtain a 
simpler expression for the functions (a, 6 = 0) 

/I sinh (X 
exp i-ax) dx = 2 , (29) 
-1 a 

which used in Eqs. (|T8|), (^, (^ and (^ leads to the parametric Ed- 

dington factor: 

/e = Jn = - - cotha, (30) 
a 

. 2/1 \ 
X = 1 coth a . 



a \a 



Eqs. (pOD were obtained previously by Minerbo M in one of the first attempts to 



obtain a variable Eddington factor for anisotropic nonequilibrium radiation from 
probabilistic arguments. Contrarily to the Lorentz case, it corresponds to a true 
nonequilibrium situation: the distribution of photons (|14D with 6 = cannot be 
transformed into a Planckian by a Lorentz transformation. 



In Figure |l] we compare the Eddington factors (pS]) , ( pOD , (p6|) and an expres- 
sion arising from a Chapmann-Enskog calculation and expressed in a parametric 
form by [0]: 

fE = 1/m — cothm, (31) 
X = — cothm ■ (1/m — cothm) 



Note that the pure heat fiux case, Eq. ( p6D grows with the energy fiux more 
rapidly than those of Lorentz and Minerbo as any advective contribution to the 
energy flux has been subtracted by demanding that /at = 0. 
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The entropy density is also modified due to the external fluxes, yielding 



s{e, Je, Jn) 



Stt 



4-^) ^^3(a,6) + ^^2(a,&) 



(32) 



{hcf3)^ 

In the limit of vanishing fluxes, a — > 0, 6 ^ we recover the equilibrium en- 
tropy multiplied by a factor QO/vr^ ~ 0.92 which is due to the approximation of 



neglecting the —1 term in the distribution function (jM]). 

In the case of pure heat flux Eq. (|32D takes a more simpler form. Using 
equations (0) and (ESf) one can express a and 6 as a function of Je and we can 



write for the entropy 



s(e, Jn 



327r 



■/cB^3(a,6) =AkBp. 



(33) 



The relation coming from the last equality is also known to hold in equilibrium. 
The nonequilibrium entropy density (|3^) in the Minerbo's limit 6 = 0, is given 



by 



(e, Jn) = ttA-Sro (e) ^ 



3/4 



27^^ '' Vsinha 
and the energy and photon densities are given by 

A8tt sinh a 



sinh a 
5 cosh a 



(34) 



{hcf /?4 a ' 
IGir sinh a 



(35) 



P 



[hcpf a ' 

Note that, the larger the flux, the lower the entropy, revealing the larger order 
existing in the system. In the limit a oo, i.e. all photons moving collectively, 
the entropy s ^ — oo and j3 oo, while in the limit a ^ we recover the 
equilibrium situation. 
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Once the nonequilibrium entropy density is known, the thermodynamic pres- 
sure can be obtained from the relation : 

s = I + I + ksl- Je + ksK ■ 4, (36) 

(being T = l/ksP) and it verifies: 

as in equihbrium. However, now the equihbrium relation between the thermody- 
namic pressure p and the trace of the pressure tensor, namely Tr(P£;) = 3p, is 
not satisfied. This can be shown by noting that the trace of the pressure tensor 
is Tr(P£;) = e for relativistic particles, and that the relation between the ther- 
modynamic pressure p and the energy density e as obtained from Eqs. (^0D, (|37|) , 



e^3(a,6) e b 
3^4[a,b) 3 3c 

differs from the equilibrium one, p = e/3. From Eq. (^), one obtains that the 
relation p = e/3 holds both in the Minerbo's case (6 = 0) and in equilibrium at 
rest {Je = 0). 

The thermodynamic pressure p is not related either to the isotropic part of 
the pressure tensor, namely e(l — x)/2, in the general situation as shown by 
inspection of Eqs.(pO|), (0) and (^): p = e(l — x)/2 only in the Lorentz's case 
and in equilibrium at rest. 

If we restrict ourselves up to second order in Je and Jn, we can find simpler 
analytical expressions for the thermodynamic quantities than the previous ones, 
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what allows a simpler physical interpretation. First, a simple relation between 



the Lagrange multipliers and the fluxes is seen to hold: 

/tv = - b, Je = 



a Ah 
3 ~ T' 



(39) 



and we can, thus, obtain an explicit flux- dependent caloric equation of state up 



to second order: 



l + -(5/|-8/^-/l, + 4/. 



with /3o : = 



247r 
{hcfe 



.1/4 



(40) 



Note that Tq = l/ksPo is the temperature corresponding to equilibrium radiation 



with the same energy density e. Thus, Eq. ( ^0]) shows that the nonequilibrium 
mean temperature defined as T = l/ksP is smaller than the corresponding local 
equilibrium temperature Tq, as expected from the general arguments of Landau 
which must hold for all possible nonequilibrium distributions of radiation. 
The entropy density of the system can also be written as the equilibrium one 



plus a quadratic correction in the fluxes: 

647r 



s(e, Je, J, 



N 



3 /3 
2 U' 



~ o ( 7/I " i^In ■ fE + In 



(41) 



and, as expected, the flux dependent contribution reduces the value of the entropy, 
in agreement to the fact that it corresponds to a more ordered physical situation. 
The particle density verifies 

IGtt 



P 



{hc(3o) 



1 2 (4'^^ '^^ 



(42) 



and the Eddington factor can be approximated by 



X 



3 + 5 



5/|-8/s-/^ + 4/^ 



(43) 
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Let us observe that both fluxes must be null in order to recover the isotropic 
Eddington factor x = 1/3, due to the fact that the fluxes /at, /e are considered 
as independent variables in the derivation. This situation is different from those 
encountered in the Eddington factor depending on two parameters introduced in 
[H, where /at and are considered as dependent variables and the Eddington 
factor depends on Je and a nonequilibrium chemical potential for the photons 
(which must be set by hand to zero at equilibrium). 

We should also notice that the correction to the isotropic Eddington factor is 
always positive (the bilinear form is always positive) up to second order in the 



fluxes. In addition, from Eq. ( |43| ) we can observe that the case of pure heat 
flux (/tv = 0) corresponds, among the situations we have considered, to the case 
in which the correction to the isotropic value 1/3 due to Je is higher: from Eq. 
( ^HD it is seen to be 4 times higher than that observed in the Lorentz situation ( 
/tv = 3/4/^;) and 5 times higher than that of Minerbo (/at = /e)- 



3 Anisotropic radiation with chemical potential 

Up to now, we have considered, as usual, that photons have a null chemical 
potential, so its number is undetermined. However, the use of a non-zero chemical 
potential for a photon gas has already been proposed in several (and physically 
distinct) situations. For example, in solid state physics, non-equilibrium but 
steady-state quasi-Fermi distributions had been used since the 1950's for the 

18 



electron gas. In this case, the study of a photon gas coupled with this electron 



gas via adsorption-emission processes leads |18| to a Bose-Einstein distribution 



2 

""^'^ = exp(/?6 + /Xp,/3)-l' ^^^^ 
where /ip^ is the chemical potential of the photons. Such a non-null chemical 
potential is related to stimulated emission of photons and other out of equilibrium 
situations. In astrophysics, massless particles like neutrinos and photons have 
been considered with nonzero chemical potential in nonequilibrium situations, 
when the number conserving Compton scattering dominates (see for example 
and references therein). 

From a technical point of view, the occupation number given in (|T^) leads 
to divergent integrated quantities out of the low-occupation number approxi- 
mation. This drawback can also be removed by introducing a nonequilibrium 
chemical potential. In addition, if we consider a photon gas with energy density 
e which transports an energy flux Je with a fixed particle flow Jjv, it seems quite 
reasonable to impose together with these constraints to the entropy, the constant 
mean number of particles p, being this latter imposed in order to be consistent 
with the fact of imposing a fixed value for the particle flow. This assumption im- 
plies, as in the examples considered above, that photons out of equilibrium may 
have a non null chemical potential. The occupation number which maximizes the 



entropy (O) under these constraints is: 



2 . ^ 

n = ^ ^ . (45) 

exp{/3pc + + I ■ pcc + K ■ c) — 1 
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Once the distribution function is known, one can obtain the thermodynamics 
of the system. As above, we restrict ourselves to the low occupation number 
limit in order to compare with the previous results. This limit corresponds to an 
ultrarelativistic classical gas (apart from the factor 2 related to the polarization of 
photons). Following the same procedure as in the previous sections, we can easily 
obtain in this approximation (assuming as in the previous section, that (3 > \I\c 
so that integrations converge and that / and K are parallel for simplicity). With 
these requirements we can find (a = \K\c and h = \I\c/P): 

P = -J^^ exp{-fiphp)'^3{a, b), (46) 

e = (^i^^^Y^ exp(-/ip,,/5)^4(a, &), (47) 

247rc 1 
Je = ^^^^^3^ exp(-/ipfe/3)- [^3ia,b) - ^iia,b)] , (48) 

Jn = j^^expi-fiphP)^ [^2(a,&) - ^3(a,&)] • (49) 
On the other hand, the entropy density can be written as 

s = kBp 



(50) 



and the thermodynamic pressure is also given by Eq. (|38[). 

With the exception of pressure, all these quantities differ from the ones in- 
troduced in the previous section. However the reduced variables, namely fE, 
/iv and X ^'^^ the same, as the new contribution due to chemical potential van- 
ishes. Therefore, this new approach does not affect the study of radiative transfer, 
though it clearly modifies the thermodynamics of the system. Note, however, that 
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the new non null nonequilibrium chemical potential for photons can be written, 
in terms of the particle density of the system as 

1 / Stt \ 

so it does not reduce to zero for null fluxes, unless we also require indepen- 
dently that p be given by the equilibrium expression p = IQn / {h(3cY . This fact 
arises directly from the hypothesis that a fixed particle density can be fixed in- 
dependently from the fluxes. However, this is also the case, for instance, in the 
treatment performed by Fu 

In order to compare with the results in the previous section, we can rewrite 
the particle density p as 



where a is the modification with respect to the particle density in (pT|) , so the 
rest of the previous quantities can also be found to be given by 

247rr 1 

Je = ^^^(1 + (^3(a, b) - ^,ia, b)) , (54) 
Jn = ^^^(1 + (^2(a, b) - ^3(a, b)) , (55) 

so that, if the chemical potential is not null, the thermodynamic equations of 
state introduced in Section 3 will be modified by a factor (1 + a). 
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4 Conclusions 

In this paper we have apphed information theory based nonequihbrium statistical 
mechanics to describe anisotropic radiation in the so-called variable Eddington 
factor scheme, investigating the main thermodynamic features of the system. 

Typically, in the variable Eddington factor closure scheme, one considers the 
reduced energy flux Je as a new variable and the pressure tensor is written as 
a function of it to close the hierarchy of radiative transfer equations. However, 
this procedure, widely used in the applications of the radiative transfer equation, 
has some thermodynamic consequences which are usually skipped. The main 
consequence is that the local equilibrium hypothesis must be abandoned as it is 
not expected to provide any dependence in the energy flux for the pressure tensor 
except in the case in which we consider radiation in a moving frame. The local 
equilibrium hypothesis implies a local black body distribution for radiation, but 



it is known [|T^ that in many radiation transfer problems the photon distribution 
strongly departs from this simpler black body expression. 

The use of information theory allows one to obtain a distribution function 
for the photon gas consistent with this closure scheme, and the fiux-dependent 
Eddington factor as well as the thermodynamics of the system can be obtained. 
In this formalism, the flux dependence of the Eddington factor merely reflects the 
fact that the fluxes characterizing the anisotropy of the system are incorporated 
as a new thermodynamical variables. The necessity of incorporating dissipative 
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fluxes within the set of thermodynamic variables beyond the local equilibrium 
hypothesis has been introduced by recent nonequilibrium thermodynamic theories 



, 21, 22 



In the variable Eddington closure scheme we have two fluxes in the system: the 
energy flux Je and the particle flux J^. If these moments are taken as dependent 
variables, we must use only one of them as a constraint in the informational 
entropy. The effect of considering each of them as a constraint has already been 
performed H, |^, |[, whereas the effect of considering them as independent 
variables is the subject of Sections 2 and 3. 

If the energy flux Je is taken as a constraint in the maximization of the infor- 
mational entropy, the results obtained by a purely macroscopic method in ||, |7| 
are recovered. As already pointed in 0, the results obtained in |^ || were 
not really a description of nonequilibrium radiation, but merely equilibrium ra- 
diation as observed from a moving reference frame. Thus, these procedures allow 
the study of anisotropic radiation but, being the anisotropy due to the relative 
motion between matter and radiation (like in the case of the cosmic background 
radiation), they do not lead to a nonequilibrium situation. Both information 
theory and the macroscopic approach followed in 0] consider an entropy de- 
pending on the energy flux but no restrictions about the global motion of the 
system are imposed. Therefore, when the condition of maximum entropy is used, 
an equilibrium moving system appears because equilibrium situations have the 
maximum entropy and the moving system verifies the imposed constraint of non 
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zero energy flux. In fact, irreversibility is related to positive entropy production 
and this fact, intrinsically related to the nature of the processes occurring in the 
system, cannot be changed by simply performing a Lorentz boost. 

Nevertheless, from the point of view of nonequilibrium thermodynamics be- 
yond local-equilibrium, we can observe that, according to the results in the 
existence of a fiux-dependent temperature might seem a physically reasonable 
assumption in this context, as it arises from a purely equilibrium situation. This 
would suggest a further validity of fiux-dependent equations of state for a more 
general, nonequilibrium situation. 

To obtain a true nonequilibrium situation, one can impose the particle flux 
Jn as a constraint, instead of the energy flux Je as done by Minerbo 0. His 
definitions were frequency-dependent and ours are not. However, if we consider 
a gray medium and the radiation intensity factorizes into a frequency-dependent 
and an angular dependent part, Minerbo's Eddington factor can be recovered 
applying information theory (with integrated frequencies) submitted to a fixed 
particle flow. 

In Section 2, by introducing both fluxes J^, as independent variables, we 
have obtained a unified treatment that allows the re-derivation of both Lorentz's 
and Minerbo's cases in the proper limits. In addition, we obtain new more gen- 
eral forms for the (two-parameter) Eddington factor that may be useful when the 
matter velocity is not null with respect to the photon flow and such an effect in- 
troduces an advective contribution to the energy flux. Hence, we have considered 
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the Eddington factor describing radiation under a pure heat flow (particle flow 
is set to zero), for which the anisotropic effects due to the flux seem to be higher 
than in any of the previous situations. 

In addition to nonequilibrium Eddington factors, we have also studied the 
nonequilibrium equations of state arising from such formalism. Note that, whereas 
the equation 

p = pksT (56) 

always remains, the nonequilibrium temperature T is related to energy by: 

1 / 247r \'/' 

so a nonequilibrium, flux- dependent temperature appears. This fact is in com- 



plete agreement with the predictions of some macroscopic theories |2^, ^] that 
have introduced nonequilibrium equations of state. Note as well that ^1^4(0, h) > 1, 
so the nonequilibrium temperature T is smaller than the equilibrium one. With 
regard to the thermodynamic pressure p, given by 

it is not related in the general case neither to the trace of the pressure tensor, 
as usual, nor to its isotropic part. This situation is, in fact, analogous to that 
encountered in for a classical ideal gas submitted to a heat flux. However, if 
b = 0, p = e/3 and it is thus related to the trace of the pressure tensor; whereas 
in the case of an equilibrium moving system (a = 0), p = e^^, i.e. it is given by 
the isotropic part of the pressure tensor. 
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Finally, we have also considered the case in which a nonequilibrium chemical 
potential for photons is introduced as suggested in [18] and done in [Q. Although 
the necessity of such an assumption is not clear, it seems a plausible ansatz which 
also removes the technical difficulties of convergence beyond the low occupation 
number limit. Within this limit, if fiph is introduced, the radiative properties of 
the system are not modified, but in the thermodynamic equations of state a new 
factor (1 + a) appears. The general case deserves further study and will be the 
object of a future paper. 
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Figure 1: Functional relationship between % and fs for four different models. 
Solid line: pure heat flux; dashed hne: Levermore (eq. (72)- (73)); short dashed 
line: Minerbo (6 = 0); dots: Lorentz (a = 0). 
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